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ABSTRACT 

We present an interesting Sunyaev-Zel'dovich (SZ) detection in the first of the Arcminute 
Microkelvin Imager (AMI) 'blind', degree-square fields to have been observed down to our 
target sensitivity of 100/xJy/beam. In follow-up deep pointed observations the SZ effect is de- 
tected with a maximum peak decrement greater than 8 x the thermal noise. No corresponding 
emission is visible in the ROSAT all-sky X-ray survey and no cluster is evident in the Palo- 
mar all-sky optical survey. Compared with existing SZ images of distant clusters, the extent 
is large (« 10') and complex; our analysis favours a model containing two clusters rather than 
a single cluster Our Bayesian analysis is currently limited to modelling each cluster with an 
ellipsoidal or spherical beta-model, which do not do justice to this decrement. Fitting an ellip- 
soid to the deeper candidate we find the following, (a) Assuming that the Evrard et al. (2002) 
approximation to Press & Schechter (1974) correctly gives the number density of clusters as a 
function of mass and redshift, then, in the search area, the formal Bayesian probability ratio of 
the AMI detection of this cluster is 7.9 x 10^:1; alternatively assuming Jenkins et al. (2001) 
as the true prior, the formal Bayesian probability ratio of detection is 2.1 x 10^:1. (b) The 
cluster mass is Mt^2()0 = 5-5^}'3X IO^^H^qAIq. (c) Abandoning a physical model with num- 
ber density prior and instead simply modelling the SZ decrement using a phenomenological 
/3 -model of temperature decrement as a function of angular distance, we find a central SZ 
temperature decrement of —295^^5 jiK - this allows for CMB primary anisotropies, receiver 
noise and radio sources. We are unsure if the cluster system we observe is a merging system 
or two separate clusters. 

Key words: cosmology: observations - cosmic microwave background - galaxies: clusters - 
Sunyaev-Zel'dovich 
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1 INTRODUCTION 

The Sunyaev-Zel'dovich (SZ) effect is the inverse-Compton scat- 
tering of cosmic microwave background (CMB) photons from the 
hot plasma within clusters of galaxies (Sunyaev & Zel'dovich 1972, 
see e.g. Birkinshaw 1999 and Carlstrom, Holder & Reese 2002 for 
reviews). The surface brightness of an SZ signal does not depend 
on the redshift z of the cluster and the integrated signal is only 
weakly dependent on z via the angular diameter distance. Hence an 
SZ-effect flux-density-limited survey can provide a complete cat- 
alogue of galaxy clusters above a limiting mass (see e.g. Bartlett 
& Silk 1994, Kneissl et al. 2001, Kosowsky 2003 and Ruhl et al. 
2004). 

Detecting and imaging the SZ effect has gradually become 
routine since it was first securely detected by Birkinshaw, Gull & 
Moffet (1981) and first imaged by Jones et al. (1993). Until re- 
cently, SZ observations have been directed almost entirely towards 
clusters selected optically or in X-ray, for example with AMI (AMI 
Consortium: Zwart et al. 2010), AMiBA (Wu et al. 2008), APEX 
(Halverson et al. 2009), CBI (Udomprasert et al. 2004), CBI-2 
(Pearson et al. 2009), OCRA (Lancaster et al. 2007), OVRO/BIMA 
(Caristrom, Joy & Grego 1996), RT (Grainge et al. 1996), SuZIE 
(Holzapfel et al. 1997), SZA (Muchovej et al. 2011) and the VSA 
(Lancaster et al. 2005). Now, however, SZ blind surveying is under- 
way, with ACT and SPT having produced initial results (Hincks et 
al. 2010, Menanteau et al. 2010, Staniszewski et al. 2009, Vander- 
linde et al. 2010 and High et al. 2010). The Arcminute Microkelvin 
Imager (AMI) is conducting a blind cluster survey at 16 GHz in 
twelve regions, each typically one deg'^, which contain no previ- 
ously recorded clusters. The AMI cluster survey focuses on depth, 
aiming to detect weak SZ-effect signals from clusters of galaxies 
with a mass above Mj.im ~ 2 x 10'^Mr;„ where Mt,200 corre- 
sponds to the total cluster mass within a spherical volume such that 
the mean interior density is 200 times the mean density of the Uni- 
verse at the current epoch. 



Table 1. AMI technical summary. Note that the brightness sensitivity is 
highly dependent on the weighting of the visibilities - in this Table we 

assume natural weighting. 





SA 


LA 


Antenna diameter 


3.7 m 


12.8 m 


Number of antennas 


10 


8 


Number of baselines 


45 


28 


Baseline length 


5-20 m 18-1 10 m 


16-GHz power primary beam FWHM 


19.6' 


5.6' 


Synthesized beam FWHM 




«30" 


Array flux-density sensitivity 


30 mjy s" 


■1/2 3mJys-V2 


Array brightness sensitivity 


4.6 mKs" 


-1/2 16mKs-i/2 


Observing frequency 


13 


.5-18.0 GHz 


Bandwidth 




3.7 GHz 


Number of channels 




6 


Channel bandwidth 




0.75 GHz 


Polarization measured 




I + Q 



The outline of this paper is as follows. In Section 2, we give a 
brief description of the instrument, observations, data reduction and 
map making techniques. Identifying cluster candidates is described 
in Section 3 - we stress that some readers will wish to jump to the 
start of Section 3 which is an important overview of the three anal- 
ysis methods and of their assumptions. We discuss how we apply 
a Bayesian analysis to the AMI data in Section 4 and present the 
results in Section 5. 

We assume a concordance ACDM cosmology, with iira = 0.3, 
Cl\ = 0.7 and Hq = 70 kms~^Mpc~^ The dimensionless Hubble 
parameter Hjq is defined as hjo = Ho/(70 kms~^Mpc~'). All coor- 
dinates are given at equinox J2000. 



2 INSTRUMENT, OBSERVATIONS, DATA REDUCTION 
AND MAPPING 

2.1 The Arcminute Microkelvin Imager (AMI) 

Sited at the Mullard Radio Astronomy Observatory, Cambridge 
19m above sea level), AMI consists of a pair of aperture- 
synthesis interferometric arrays optimised for SZ-effect imaging 
centred at 16 GHz, with six frequency channels. The Large Array 
(LA) has a high resolution and flux-density sensitivity and is used 
primarily to detect contaminating sources which can then be sub- 
tracted from the Small Array (S A) maps. AMI is described in detail 
in AMI Consortium: Zwart et al. (2008) and the technical aspects 
of the arrays are summarised in Tabic 1 . The S A has been operating 
since 2005 (see e.g. AMI Consortium: Barker et al. 2006, Scaife et 
al. 2008 and Scaife et al. 2009) and the LA since 2008 (see e.g. 
AMI Consortium: Hurley- Walker et al. 2009). Pointed SZ observa- 
tions have been straightforward but for blind observations we have 
felt it essential to get the best control of systematics that we can - 
for example, we found hard-to-unravel problems with LA pointing 
and errors in the electrical lengths of the lags in both the LA and the 
SA Fourier-transform correlators that produce small position shifts 
- we now have corrections for these problems that are adequate. We 
also have very good control over the influence of radio source con- 
tamination (see e.g. AMI Consortium: Franzen et al. 2009, Feroz et 
al. 2009, Waldram et al. 2010 and AMI Consortium: Davies et al. 
2011). 
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2.2 Observations 

The results presented here are from observations of field AMI002 
which is centred on 02^ 59™ 30' +26° 16' 30". AMI002 is the first 
field to have been analysed as it was the first to reach a target depth 
of 100;UJy/beam. SA observations of AMI002 began on 2008 July 
19 and ran until 2010 March 3, by which time 1200 hours of data 
had been gathered; LA observations began on 2008 August 8 and 
ran until 2010 January 10, collecting 630 hours of data. Using both 
the SA and the LA the field was typically observed for 8 hours in 
a day; this often comprised two individual observations each of 4 
hours, split up with an observation of a flux-density calibrator. Ob- 
servations were started at different positions in the field to improve 
the Mv coverage. 

A rastering technique was used for both the LA and the 
SA survey observations, where the pointing centres lie on a 2-D 
hexagonally-gridded lattice. The LA observations form a pai't of the 
IOC survey data, which are described in detail in AMI Consortium: 
Franzen et al. (201 1). Additional dedicated pointings towards the 
cluster candidates are included to ensure that maximum sensitiv- 
ity was obtained in the LA maps. For the IOC survey observations, 
the pointing centres are separated by 4 arcmin, which allows us to 
obtain close to uniform sensitivity over the field while minimising 
the observing time lost to slewing. In order to detect all important 
sources within the SA field, the LA field is slightly larger and the 
thermal noise is typically a factor of two lower than the SA ther- 
mal noise. To account for the SA map noise (cTgA.survey) increasing 
towards the edge of the field, the LA map consists of two distinct 
regions, the inner and the outer. The iimer area of the LA field was 
observed to a noise level of ~ 50/i Jy, whereas the noise in the outer 
area was approximately twice as high. The outer region of the LA 
map is also used to detect bright sources lying just outside the SA 
field. The resulting LA noise map is shown in Figure 1. For the 
SA survey observations the pointing centres are separated by 13 ar- 
cmin giving a closc-to-uniform noise level of ~ lOO/iIy over the 
map. The SA noise map is shown in Figure 2. Follow-up SA obser- 
vations towards the cluster consisted of 50 hours of data centred at 
03'' 00" 08.66* +26° 15' 16.1" resulting in a noise level of 65fiJy. 

The phase calibrator was observed for two minutes every hour 
using the SA and for two minutes every ten minutes using the LA. 
The phase calibrator used for both the LA and the SA observa- 
tions was J0237+2848. The amplitude calibration for the SA uses 
3C286 and 3C48 which are observed daily; the assumed flux den- 
sities are shown in Table 2 (AMI Consortium: Franzen et al. 201 1). 
The LA was flux-density-calibrated from the SA measurements 
of J0237+2848; we have adopted this approach to minimise inter- 
array calibration errors. Although the flux density of J0237+2848 is 
known to vary at AMI frequencies, with a mean variability index of 
3.1 over 269 days (AMI Consortium: Franzen et al. 2009), we en- 
sured that SA measurements of this source were always within 30 
days of the LA observations. This calibration scheme is described 
in detail in AMI Consortium: Franzen et al. (2011) and provides 
flux-density calibration errors of less than 5%. 



2.3 Data Reduction 

There are 65 LA observations and 337 SA observations of AMI002, 
each being passed through reduce, the in-house software devel- 
oped for the VSA and AMI data reduction. REDUCE was used to 
flag telescope pointing errors, shadowing effects and hardware er- 
rors. The data are also flagged for interference before being Fourier 
transformed into the frequency domain, where they are corrected 




3'03" a'oO" 2''58' 

Right Ascension (J2000) 



Figure 1. Noise map for the LA survey field. The inner region noise is « 
50;Ujy, while the noise on the outer region is ^ lOO^lJy. The hexagonal 
region around 03^ 00™ 10' +26° 15' 00" is next to the cluster and was 
observed to f» 30/iJy. The inner region of the noise map consists of three 
subregions; these have slightly different sensitivities due to varying weather 
conditions and slight differences in observing time. 



Table 2. Assumed flux densities for the SA flux-density calibrators. 



Channel 


v/GHz 


Si+Q/Jy 






3C48 


3C286 


1 


14.2 


1.850 


3.663 


2 


15.0 


1.749 


3.535 


3 


15.7 


1.658 


3.414 


4 


16.4 


1.575 


3.308 


5 


17.1 


1.500 


3.206 


6 


17.9 


1.431 


3.111 



for system-temperature variations, phase-calibrated and amplitude- 
calibrated. In the frequency domain the data are again searched for 
interference and baselines with inconsistent flux-density values are 
flagged. The data are reweighted so that baselines and channels 
with the lowest noise have the highest weight. The data are then 
stored as uvfits files, with each raster pointing being treated as 
an independent source within the FITS definition. This reduction 
scheme follows that of AMI Consortium: Davies et al. (2009). In- 
dividual UVFITS files for the LA and the SA are combined into a 
single multisource UVFITS file for each array, which are taken into 
AiPS ' for imaging. 

The SA data were checked for systematics using two jack- 
knife tests. In test (a), calibrated visibilities from "plus" correlator 
boards are subtracted from those obtained from "minus" correla- 
tor boards - the signal is the same for both correlations but the 
latter inserts an additional 180° phase shift into the signal from 
one antenna (see Holler et al. 2007 for a full description of the 
AMI correlator). For test (b), data obtained before the weighted 
median date of the visibilities are subtracted from data obtained 

' http://www.aips.nrao.edu 
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Figure 2. Noise map for the SA AMI0()2 field. Tlie noise at tlie edge of 
the map increases due to the primary beam of the SA. In the central region 
the map noise is lOO^Jy. This noise map does not include the follow-up 
SA observations towards 03^ 00™ 08.66" +26° 15' 16.1" as those data were 
analysed separately. The noise level of the foUow-up SA observations is 
65/JJy. 

afterwards. Either test will remove signals present in both halves 
of the data but noise or systematics that vary with time will re- 
main. For the follow-up pointed SA observations presented in this 
paper test (a) revealed no systematics and test (b) showed a neg- 
ative feature with a flux-density of 0.35mJy/beam associated with 
the 2.26mJy/beam source at 03:00:29.46 +26:18:39.9. Investiga- 
tion demonstrated that this residual was a consequence of the flux- 
density of the source being dependent upon the elongation of the 
synthesized beam. For reasons of scheduling, it became clear that 
the synthesized beam from the first half of data was elongated in 
the NW-SE direction and the source was measured to have a flux- 
density of 2.5 ImJy/beam. In the second half of the data the synthe- 
sized beam was extended in the NE-SW direction and the measured 
flux-density was 2.03mJy/beam. Maps of the jack-knifed pointed 
SA observations are shown in Figure 3. 

2.4 LA map-making and source-finding 

LA maps for each AMI channel and the continuum were pro- 
duced for each of the pointings within the AMI002 field using the 
AIPS task IMAGR. The maps are CLEANed to three times the map 
thermal-noise without any individual CLEAN boxes. The individ- 
ual pointings are combined using the flatn task, discarding data 
lying outside the 10% point of the power primary beam, flatn 
is also used to create appropriately-weighted noise maps using the 
thermal noise levels in the individual pointings. 

Source finding is carried out using the LA continuum map 
with the AMI SOURCEHND software. All pixels on the map with 
a flux density greater than 0.6 x 4 x <t„, where (T„ is the noise 
map value for that pixel, are identified as peaks. The flux densi- 
ties and positions of the peaks are determined using a tabulated 
Gaussian sine degridding function to interpolate between the pix- 
els. Only peaks where the interpolated flux density is greater than 



4 X CT„ are identified as sources. The AIPS routine JMFIT fits a two- 
dimensional Gaussian to each source to give the angular size and 
the integrated flux density for the source. These fitted values are 
compared to the point-source response function of the telescope to 
determine whether the source is extended on the LA map. The map- 
ping and the source finding techniques are described in more detail 
in AMI Consortium: Franzen et al. (2011). 

For each source we use the SOURCEhnd algorithm to find 
the flux densities in the individual AMI LA channel maps at the 
positions of the detected sources. By assuming a power-law rela- 
tionship between flux density and frequency (5 0= v~") we use the 
channel flux-densities to determine the spectral index a for each 
source. The spectral index is calculated using an MCMC method 
based on that of Hobson & Baldwin (2004) - the prior on the spec- 
tral index has a Gaussian distribution with a mean of 0.5 and a of 
2.0, truncated at ±5.0. The minimum spectral index of a source in 
the AMI002 field was found to be 0.0 and the maximvun was 1.8. 
The map noise in each channel map at the position of the source 
was used to calculate the weighted mean of the channel frequencies 
and determine the effective central frequency Vq of the source. The 
effective central frequency varies between pointings due to flag- 
ging applied in reduce. Unlike in AMI Consortium: Franzen et 
al. (2011), the data are not reweighted to the same frequency be- 
cause this leads to a small loss of sensitivity. 

In total we detect 203 sources in the AMI002 LA map at four 
times the LA map noise (cfla), 1 1 of which are extended. The most 
extended source has an area of 1.9 LA synthesized beams. As the 
SA synthesized beam is significantly larger we do not expect any 
extended sources in the SA map. For each source we catalogue the 
right-ascension Xs, declination ys, flux density at the central fre- 
quency 5o, spectral index and the central frequency. If a source is 
extended we use the centroid of the fitted Gaussian as the position 
and the integrated flux density instead of the peak flux density. 



3 roENTIFYING AND MODELLING CLUSTER 
CANDIDATES 

Our analysis necessarily depends in part on the fact that we do not 
know - in the absence of e.g. optical spectroscopic observations 
- the redshift of the blind SZ clusters. We have thus carried out 
our analysis in two main ways, both fully Bayesian and based on 
Hobson & Maisinger 2002, Marshall et al. 2003 and Feroz et al. 
2009, as follows. 

(1) Physical model. We assume an isothermal jS -profile for the 
gas density as a function of radius; we assume all the cluster kinetic 
energy is in the internal energy of the cluster gas and that the rela- 
tion between gas temperature and total cluster mass is then given 
by the virial theorem; and we assume the prior probability for the 
comoving number density of clusters as a function of total mass 
and redshift is given by previous theoretical/simulation work - we 
here use the predictions of Evrard et al. (2002) and Jenkins et al. 
(2001) and note that more recent such work does not make a sub- 
stantial difference for our purposes. With these assumptions we are 
then able to (a) estimate the significance of an SZ detection, and 
(b) produce probability distributions of physical cluster paraineters 
such as mass and radius. For both (a) and (b) the methodology takes 
into account radio sources, receiver noise, and the statistical prop- 
erties of the primordial CMB structure; it cannot take into account 
other effects that have not been dealt with in, for example, telescope 
design, telescope commissioning, observing and data reduction. 

(2) Phenomenological model. Some or all of the assumptions 
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Figure 3. The jack-knifed SA pointed observations of the cluster candidate (an image of the data is shown on the left of Figure 6). On the left, data are spUt 
into plus and minus baseline, one the right, data are spUt according to median date. The contour levels are linear from 2 asA,sttrvey to 10 OsA,survey (<^SA,pointed 
= 65/iJy); positive contours are solid lines and negative contours are dashed Unes. The ellipse at the bottom left of the maps shows the SA synthesised beam. 



in (1) may be poor or wholly wrong. Accordingly in (2), we make 
far fewer asswnptions. We assume isothermality and that the tem- 
perature decrement as a function of angular distance is given by a 
/3-modeI. This model cannot give probability distributions of values 
of physical importance such as mass, but still does give the signifi- 
cance of the SZ detection in the presence of radio sources, receiver 
noise and primordial CMB structure; like (1) it cannot take into 
account other effects that have not previously been dealt with. 

We give the significance of decrement detection in a third way, 
the decrement signal in units of receiver noise. We point out that for 
AMI this method of course takes no account of primordial CMB 
structures but does take into account radio sources, and the higher 
flux-density sources have had their SA flux-densities estimated in 
a Bayesian way from SA data and priors from LA measurements 
- this allows for inter-array calibration problems and for LA and 
SA observations that were not precisely simultaneous. Again, Uke 
(1) and (2), this method cannot take into account other effects that 
have not previously been dealt with. 



3.1 Physical model 

Our primary Bayesian analysis is based on a physical model 
for the cluster producing the SZ effect. The SA observa- 
tions of the AMI002 survey field are analysed using a model 
characterised by the parameters = (©c,*?), where = 
{xc,yc,<j>,f,P,rc,Mj20QjZ) are cluster parameters and I* — 
{xs,ys,So,a) are source parameters (Feroz et al. 2009). Here Xc 
and yc give the cluster position, is the orientation angle measured 
from N through E, / is the ratio of the lengths of the semi-minor to 
semi-major axes, P describes the cluster gas density Pg according 
to Cavaliere & Fusco-Femiano (1976,1978), where the gas density 
decreases with radius r 



Pgir 



Pg(0) 



[l + (r/r 



(1) 



which the mean total density is 200 times the critical density Pcnt. 
Feroz et al. (2009) and AMI Consortium: Rodriguez-Gonzalvez et 
al. (201 1) describe the parameters and the methods used to extract 
these from the data in more detail. For this work we sample from 
Xc, yc, (l>,f,P, rc, Mj- 200 and z and derive other cluster parameters 
such as the cluster gas mass Mg 2oo aid the cluster temperature T. 
We also assume a mass-temperature relationship characteristic of a 
virialised cluster; this is the favoured model (M3) in AMI Consor- 
tium: Rodriguez-Gonzalvez et al. (201 1), although we sample from 
Mt.ioq rather than Mg 200, see also AMI Consortium: Olamaie et al. 
2010. The total cluster mass within r2oo is 

Mt,200 = y r^oo(200pcrjt). (2) 

The gas fraction fg is derived from the results of Komatsu et al. 
(2010) taking into account our value for h and that the gas-mass 
fraction is ~ 0.9 of the baryonic mass fraction. The ellipticity of 
the clusters is calculated by applying a coordinate transformation 
from point {61 , G2) on the sky: 

l/\/7 J \ -sin0 cos(j> )\ 02 



(3) 



rc is the core radius, Mj 200 is the cluster total mass within a radius 
r2oo and z is the cluster redshift. r2oo is defined as the radius inside 



Lines of constant 6' represent ellipses enclosing an area Kab, where 
a is the semi-minor axis and b is the semi-major axis. This trans- 
forms the circular slices perpendicular to the line of sight to an el- 
lipse, keeping the area of the ellipse the same as the circular slice. 
A summary of the priors used on the model parameters is shown in 
Table 3. 

The above approach has already been used to detect the SZ ef- 
fect from AMI observations of known clusters in AMI Consortium: 
Zwart et al. (2010) and AMI Consortium: Rodriguez-Gonzalvez et 
al. (201 1). However, for blind cluster surveys we are faced with the 
additional problem that we do not have a priori evidence for a clus- 
ter at a particular position (or redshift). In analysing a survey field, 
the marginalised posterior distribution in the (xe.je) -plane will typ- 
ically contain a number of local peaks; some of these may corre- 
spond to the presence of a real cluster, whereas others may result 
from chance statistical fluctuations in the primordial CMB and/or 
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Table 3. Priors used for the Bayesian analysis assuming a physical cluster model. 



Parameter 



Prior 



Source position (Xj) 
Source flux density {So /Jy) 
Source spectral index (a) 
Redshift (z) 

Core radius (r^/A^Q^kpc) 
Beta (i3) 

Mass {MT,20o/h^oM@) 
Gas fraction (Jg/h^Q) 
Cluster Position (x^) 
Orientation angle (0 / deg) 

Ratio of the length of semi-minor to semi-major axes (/) 



A delta-function prior using the LA positions 
A Gaussian centred on the LA continuum value with a <7 of 40% 
A Gaussian centred on the value calculated from the LA channel maps with the LA error as (7 
Joint prior with Mt between 0.2 and 2.0 (Jenkins et al. 2001 or Evrard et al. 2002) 
Uniform between 10 and 1000 
Uniform between 0.3 and 2.5 
Joint prior with z between 2.0 x lO"* and 5 x lO" (Jenkins et al. 2001 or Evrard et al. 2002) 
Delta-function prior at 0.1 1 (Komatsu et al. 2010) 
Uniform search triangle (Figure 4) 
Uniform between and 180 
Uniform between 0.5 and 1.0 



instrument noise. Each local peak in the posterior is automatically 
identified by the MultiNest sampler (Feroz & Hobson 2008 and 
Feroz, Hobson & Bridges 2008) used in our Bayesian analysis, and 
may subsequently be analysed independently to obtain cluster pa- 
rameter estimates. 

To determine the significance of each such putative cluster de- 
tection, we perform a Bayesian model selection, which makes use 
of estimated cluster number counts from analytical theory (e.g. the 
Evrard et al. 2002 approximation to Press & Schechter 1974) and 
numerical modelling (e.g. Jenkins et al. 2001) together with mea- 
surements of the rms mass fluctuation amplitude on scales of size 
8 AjQ^Mpc at the current epoch (see e.g. Lahav et al. 2002, Seljak 
et al. 2005 and Vikhlinin et al. 2009). It must be borne in mind, 
however, that the actual values of the number density of clusters, 
particularly at high redshift, are uncertain and hence the degree of 
applicability of these as priors is unclear. 

In our Bayesian model selection, we calculate the formal 
Bayesian probability of two hypotheses: the first, Pr(//>i|D), as- 
sumes at least one cluster with Mj ^aa < ^r,200 < ^r.max is asso- 
ciated with the local peak in the posterior distribution under consid- 
eration; the second, Vx{Hq\D), assumes no such cluster is present. 
Here Mj \^^ is the limiting cluster mass that can be detected and 
M-T^max is the maximum mass of a cluster. In particular, we con- 
sider the ratio R (also known as the Bayes factor, or the odds) of 
these two formal probabiUties 



R 



_ Pr(//>i|D) 
Pr(/fo|0) ■ 



(4) 



To evaluate this ratio, let us first denote by 5 the area in the {xc,yc)- 
plane of the 'footprint' of the local posterior peak under consider- 
ation (we will see below that a precise value for S is not required). 
Also, we denote by H„ the hypothesis that there are n clusters with 
^r,Um < Mj 2m < Mj^max with centres lying in the footprint 5, so 
that 



n=l 

Thus equation (4) can be written as 



ir=lPr(g«|0) _ L:=iPr(D\H„)Pr(H„) 



Pt{Hq\D) 



Pr(D|tfo)Pr(Ho) 



(5) 



(6) 



where we have used Bayes' theorem in the second equality. Assum- 
ing that objects are randomly distributed over the sky, then 



Pr(H„) 



n! 



(7) 



where Hs is the expected number of clusters with Mrum < 
Mt.ioq < Afr.max in a region 5. This is given by Hs = 5jU, where 
H is the expected number of clusters per unit sky area: 



M= / / 

-'Zmn JMr 



V.max tfin 



dMdz 



dMdz, 



(8) 



where n{z,M) is the comoving number density of clusters as a func- 
tion of redshift and mass. For the calculation of /x, we follow the 
method of either Evrard et al. (2002) or Jenkins et al. (2001). If we 
further assume that there is very low probability of two or more 
clusters having their centres in the region S (jls <C 1) we can ne- 
glect ju| and larger powers of Hs, so that equation (6) can be ap- 
proximated simply by 



R! 



Zo 



(9) 



where the Zi(5) = Pr(£)j/fi) is the 'local evidence' (see Feroz et 
al. 2009) associated with the posterior peak under consideration in 
the single-cluster model, and Zq = Pr(D|i/o) is the 'nuU'evidence 
(which does not depend on 5). 

Our Bayesian analysis uses MultiNest to calculate the 
Bayesian evidence for the different hypotheses (Feroz & Hobson 
2008 and Feroz, Hobson & Bridges 2008). When searching for 
clusters in some survey areaA, however, a uniform prior n{xc,yc) = 
1 /A is assumed on the position of any cluster, rather than assuming 
a uniform prior over the footprint S. Thus, MultiNest returns a 
local evidence associated with the posterior peak that is given by 



Zi(5) = -Zi(5), 



(10) 



and the 'null' evidence Zq = Zq remains unchanged. Thus, if we 
denote the expected number of clusters in the survey area by jUa = 
{A/S)ns, then (6) becomes 



R'. 



(11) 



Here Z\ (5) and Zq are outputs of MultiNest and Ha is easily 
calculated from (8) given some assumed cluster mass function, and 
so R may then be calculated. In our analysis /X/i < 1 and the R 
value that we calculate is smaller than that obtained by setting the 
prior ratio equal to unity. Jeffreys (1961) provides an interpretive 
scale for the R value, as do revised scales such as Gordon & Trotta 
(2007). Moreover, the R value in (11) can be turned into a formal 
Bayesian probability p that the putative detection is indeed due to 
a cluster with mass Mj^m < Mt;ioo < ^r,max and centre lying in 
S, which is given by 
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Table 4. Priors used for the Bayesian analysis assuming a phenomenologi- 
cal cluster model. 



Parameter 


Prior 


Aro 


Uniform between ±3000^K 


ec 


Uniform between 20" and 500" 


P 


Uniform between 0.4 and 2.5 


R 




^ l+R' 





(12) 



3.2 Phenomenological model 

An alternative approach is to set aside the physical cluster model 
and instead adopt a model based on a phenomenological descrip- 
tion of the SZ decrement itself. 

In this case, at the location of each putative cluster detection 
identified using the physical cluster model, we simply fit a j3 pro- 
file to the SZ temperature decrement using the parameters dc, /3 
and ATq to characterise shape and magnitude of the decrement ac- 
cording to 

ATsz = AToi^l + -^j . (13) 

The assumed priors on these parameters are summarised in Table 4. 
In this analysis we continue to use Gaussian priors on the flux den- 
sities and on the spectral indices of significant sources, and delta- 
function priors for faint sources. We also assume a Gaussian prior 
(<T = 60") on position centred on each decrement. 

This approach allows us to produce a posterior distribution 
that directly describes the temperature decrement and also allows 
us to evaluate what proportion of the decrement is caused by the SZ 
effect, while also accurately accounting for point sources, receiver 
noise and the statistical properties of the primary CMB anisotropy. 




3'"oi" a'oo" s'so" 

Right Ascension (J2000) 



Figure 4. Noise map for a SA triangle of observations out to the 0. 1 con- 
tour of the power primary beam. The inner triangle is between the pointing 
centres; the outer triangle is the area that is searched for clusters with our 
Bayesian analysis. 



We follow up our most significant detections with pointed ob- 
servations towards the candidate. The data from these observations 
can be analysed with our Bayesian method with lower dimensional- 
ity because there are fewer sources within 0. 1 of the power primary 
beam with flux densities greater than 4(TsA,pointed- For the follow- 
up pointed observations the prior on the cluster position is altered 
to a 1000" X 1000" box centred on the pointing centre and we al- 
low our Bayesian analysis software to fit the source positions with 
a Gaussian prior centred on the LA position with an error of 5". 



4 THE ANALYSIS 

The AMI002 SA survey map contains 24 individual pointing cen- 
tres. A single Bayesian analysis of the entire field is prohibitively 
computationally expensive because of the large quantity of data and 
the high dimensionality of the parameter space. Instead three point- 
ings are analysed at a time. Each set of three pointing centres form 
a triangle and in total there are 30 different triangles in the AMI002 
field, an example of which is shown in Figure 4. 

To reduce the dimensionality of the parameter space further, 
all sources located at positions where the primary beam has fallen 
below 10% of its maximum, together with sources that have a flux 
density measured on the LA that is lower than 4crsA.suivey! ars given 
delta-function priors on their positions, spectral indices and flux 
densities. We search for clusters in a triangular area which is an en- 
larged version of the triangle formed between the pointing centres - 
the radius of the inscribed circle is 3' larger This allows us to detect 
clusters out to the edge of our most sensitive areas and ensures that 
the search areas belonging to adjacent triangles overlap. The min- 
imum rms noise within a search triangle in the AMI002 field is ~ 
100/iJy and the maximum is ~ 140/xJy. The limiting cluster total 
mass is set to Mjun, = 2 x lO'^h^^'Mp) and the maximum cluster 
mass to Mj jnax = 5 X lO'^hyjj Mpj. The limiting mass is conserva- 
tive given the radio flux-density sensitivity of our observations. 



5 RESULTS AND DISCUSSION 

The most significant candidate cluster detection made using the 
Bayesian analysis of the AMI SA survey field AMI002 is located 
at I 03*' 00™ 16.5' +26° 13' 59.5", where MultiNest identifies a 
single marginalised posterior peak in the (.TcO'c) -plane centred on 
this location. The significance of the cluster detection is /? ~ 8.7 
when we use Model (1) and the Evrard et al. (2002) prior and 
~ 26 when we apply Model (1) using the Jenkins et al. (2001) 
prior The relevant area of the survey field is shown in Figure 5 be- 
fore and after source subtraction (see below). In the search triangle 
that contains our cluster candidate there are 59 sources within 0.1 
of the power primary beam, 43 of which have a flux density below 
4fsA.survey; the Other 16 have been modelled with our Bayesian 
analysis. The location of the marginalised posterior peak is indi- 
cated by the small box in the figures. 

At this position in our survey field is a highly-extended, non- 
circular negative feature with a peak flux-density decrement of ~ 
0.6 mJy (5cTsA. survey)- SA observations are mapped in AIPS using 
the same method as for the LA, but with a pixel size of 15arcsec. 
We subtract sources from the mv-FITS data using the in-house soft- 
ware MtJESLI. MUESLI performs the same function as the AIPS task 
UVSUB; however, it is optimised for processing AMI data. The pa- 
rameters of the 16 modelled sources are shown in Table 5; we find 
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Figure 5. The SA survey-field observations of the cluster candidate. On the left is the map before source subtraction and on the right is after the map after the 
sources in Table 5 have been subtracted. The + symbols indicate the positions of sources with flux densities less than 4 C7sA,survey and the x symbols represent 
sources which have a flux density greater then 4 OsA.survey The box on the source-subtracted map shows the position of the cluster candidate. These maps 
are signal-to-noise maps. Each pointing has a thermal-noise level of approximately 0.11 mjy. The contour levels are linear from 2 asA,survey to 10 (ysA,survey\ 
positive contours are soUd lines and negative contours are dashed lines. The ellipse at the bottom left of the maps shows the SA synthesised beam. 



no evidence that any of them is extended relative to the LA syn- 
thesised beam. The source subtraction leaves very little residual 
flux density on the map, indicating that the phase stability and cali- 
bration of AMI is robust. The most significant source-subtraction 
residuals are towards the edge of the SA power primary beam 
where we expect the phase errors to be larger and the beam model 
to be less accurate. 

The cluster candidate was followed up with a pointed obser- 
vation. Within the 10% point of the SA power primary beam 31 
sources were observed with the LA, 9 of which were detected at 
above 4 CgA.pointed ™d are modelled by our Bayesian analysis. 
These 9 sources are a subset of the 16 sources modelled on the 
Bayesian analysis of the survey data; they are indicated by a 'tick' 
in the last colimm of Table 5. We find no evidence that any of these 
9 sources is extended relative to the LA synthesised beam. The im- 
age produced from the pointed-observation data is shown before 
and after source subtraction in Figure 6. Again we see a highly- 
extended, non-circular negative feature with a peak flux-density 
decrement of ~ 0.6 mJy (ScTsA.pointed)- The SA synthesized beam 
and the uv coverage for the follow-up pointed observation is shown 
in Figure 7. To estimate the maximum level of contamination from 
the residuals of the sources in Table 5, we assume that the residual 
is equal to the error in the source flux and sum the absolute value of 
the synthesized beam contribution from each of these residuals at 
the positions of candidate 1 and candidate 2, we find contributions 
of 32;Ujy and 70/iJy respectively. Hence, if in the unlikely case 
all sources leave a feature of magnitude equal to the error in that 
source flux, and that these features conspire in such a way to con- 
tribute only negative flux at the positions of candidates 1 and 2, we 
find the contribution to the total SZ signal is minimal. This calcula- 
tion does not account for any errors in the shape of the synthesized 
beam, due to e.g. antenna positions. 

Our Bayesian analysis of the pointed-observation data, which 
have a higher signal-to-noise ratio than the survey data, finds two 
local peaks in the marginalised posterior distribution in the {xc,yc)- 



plane. These cluster candidates are: candidate 1 at J 03*' 00" 14.8' 
-1-26° 10' 02.6" and candidate 2 at J 03'> 00"" 08.9' +26° 16' 29.1". 
The Model (1) significance of the two cluster detections are R\ ~ 
7.9 X lO'* and R2 ~ 560 respectively when we apply the Evrard et al. 
(2002) model and Ry^lAx 10^ and R2 ~ 1800 respectively when 
we apply the Jenkins et al. (2001) model. These larger values for 
the /?-ratio, as compared with those obtained using the survey data, 
result from the higher signal-to-noise ratio of the pointed observa- 
tion. The evidence values, 7J-ratios and related parameters for the 
survey observations and the pointed observation are summarised in 
Table 6. We also made a direct comparison of the Bayesian evi- 
dence for a model containing two clusters and a model containing 
just a single cluster and find that the Bayesian evidence is 7.6 x 
10^ higher for the model containing two clusters. 

The ID and 2D marginal posterior distributions for a selection 
of the physical parameters of each cluster are shown in Figure 8. 
We are able to constrain 200. even though it is dengenerate with 
z. As Mr 200 is large this degeneracy causes the derived z value to be 
low. We are able to constrain the well known degeneracy between 
jS and rc and find that values of /3 < 1.0 do not fit our data. We 
also find that the best-fit ratio of the lengths of the semi-minor to 
semi-major axes is 0.6 and 0.75 respectively; the orientation angles 
are 122° and 78°. 

Finally, we used Model (2) and performed a Bayesian analy- 
sis where we simply fit a j3 profile to the SZ temperature decrement 
directly. The ID and 2D marginal posterior distributions for the pa- 
rameters 6c, jS and ATq are shown in Figure 9. From our data we 
are able to tightly constrain ATq at ^ — 300/xK, but we are unable 
to accurately derive j8. The mean values and 68% confidence Umits 
for each parameter are given in Table 7 and demonstrate signifi- 
cance of the detections directly. 

We have looked for optical identification of the cluster in the 
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Table 5. The fitted parameters for the 16 sources with LA flux densities greater than 4 ffsA.survey (65jU/y). This includes the nine sources with flux densities 
greater than 4(TsA,pomted- The positions and mean frequencies are from LA observations, whereas the flux densities and spectral indices are the values obtained 
from our Bayesian analysis of the SA survey field. 



Right ascension 


Declination 


Flux density 


Spectral index 


Mean frequency 


Flux density > 4asA,poiiited 


(J2000) 


(J2000) 


(mjy) 




(GHz) 






4-96-1 Q-4n 8'^ 


1 21 ± 12 


+ i.jo HI U.JO 


15 63 


/ 

V 


03:00:29.46 


+26:18:39.95 


2.26 ±0.12 


+0.71 ± 0.29 


15.64 


V 


02:59:06.92 


+26:15:29.59 


0.26 ± 0.09 


+1.46 ± 1.10 


15.57 


V 


02:59:50.35 


+26:25:22.37 


0.23 ±0.11 


+0.51 ± 1.24 


15.58 


X 


02:59:39.76 


+26:05:56.15 


0.40 ± 0.09 


+2.31 ± 1.06 


15.52 


X 


03:00:15.23 


+26:19:25.56 


1.44 ±0.11 


+1.59 ± 0.40 


15.64 


V 


02:59:23.57 


+26:05:54.53 


0.40 ±0.10 


+0.81 ± 1.24 


15.54 


X 


02:59:55.16 


+26:27:26.24 


8.49 ± 0.22 


+0.33 ± 0.07 


15.59 




02:59:10.71 


+25:54:31.60 


3.81 ± 0.41 


+1.05 ± 0.16 


15.57 


X 


02:59:35.43 


+26:17:26.77 


0.64 ± 0.09 


-0.22 ± 0.94 


15.53 


V 


03:00:49.28 


+26:15:05.70 


0.53 ± 0.12 


+0.40 ± 0.42 


15.67 




02:59:29.68 


+26:09:46.99 


0.50 ± 0.07 


+1.71 ± 1.26 


15.55 


X 


02:59:41.05 


+26:02:20.41 


1.58 ±0.12 


+1.36 ±0.36 


15.54 




02:59:57.17 


+25:53:56.17 


1.40 ± 0.23 


+2.20 ± 0.34 


15.58 


X 


03:00:01.33 


+26:21:01.55 


1.96 ±0.12 


-0.45 ± 0.31 


15.56 


V 


02:58:25.32 


+26:16:59.59 


1.63 ± 0.33 


+0.95 ± 0.27 


15.58 


X 




03 01 00 00 45 30 15 0002 59 45 30 15 03 01 00 00 45 30 15 0002 59 45 30 15 

RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

Cont peak flux = 2.681 1 E-03 JY/BEAM Cont peak flux = -6.1271E-04 JY/BEAM 

Leva = 6.650E-05 * (-10, -9, -8, -7, -6, -5, -4, Leva = 6.650E-05 * (-10, -9, -8, -7, -6, -5, -4, 

-3, -2, 2, 3, 4, 5, 6, 7, 8, 9, 10) -3, -2, 2, 3, 4, 5, 6, 7, 8, 9, 10) 



Figure 6. The AMI SA pointed observation towards the cluster. On the left is the map before source-subtraction and on the right is the map after the relevant 
sources in Table 5 have been subtracted. The + symbols indicate the positions of sources with flux densities less than 4 asA,pomted> the x symbols represent 
sources which have a flux density greater than 4 C7sA,pomted on the SA map. For this run, the sampler has been allowed to fit the positions of the x type sources 
with a Gaussian centred on the LA source position. The contour levels are linear from 2 asA.pointed to 10 CJsA.pointed (OsA,pomted = 65|tJy); positive contours are 
soUd lines and negative contours are dashed lines. The boxes indicate the positions of the cluster candidates. Candidate 1 is at J 03* 00" 14.8* +26° 10' 02.6" 
and candidate 2 is at J 03'^ 00" 08.9^* +26° 16' 29.1". When imaging the source subtracted map clean boxes have been placed around each candidate. The 
eUipse at the bottom left of the maps shows the SA synthesised beam. 



Table 7. Mean values and 68% confidence limits for the parameters in the 
SZ decrement model for candidate 1 and candidate 2. 



Parameter 


Pointed (candidate 1) 


Pointed (candidate 2) 


ed" 
P 

ATo/jUK 


i56;g 

-295+lt 


121+'' 
1 46+ '"3 

-302+™ 



Palomar all-sky survey and X-ray identification from ROSAT ^ - 



no cluster identification is evident. We plan to perform X-ray and 
optical follow-up observations. 



6 CONCLUSIONS 

• We have presented a large, complex Stmyaev-Zel'dovich 
structure in an AMI blind field. The structure may be two separate 
components or be a single merging system. 

• A Bayesian analysis using a physical model for the cluster (in- 
cluding assimied priors on the number density of clusters) was used 
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-3, -2, 2, a, 4, 5, 6, 7, 8, 9, 10) Distance (arcmin) Kilo Wavlngth 



Figure 7. On the left is the synthesized beam for the AMI SA pointed observation towards the cluster (Figure 6). The contour levels range from 6% to 30% 
with intervals of 3%; positive contours are solid lines and negative contours are dashed lines. The image in the centre shows one-dimensional slices through 
the centre of the synthesised beam; dashed Hnes show the profile along the declination axis and solid lines show the profile along the right ascension axis. The 
image on the right shows the con'esponding uv coverage of the observation; a different colour is used for each AMI channel. 



Table 6. Evidences, R-ratios and related parameters for the detection of the cluster candidates in the triangle of survey observations and the follow-up deep 
pointed observation. The limiting total mass is Mj- ii,,, = 2 x lO''* h^^Mcj. 



Parameter 


Survey 


Pointed (candidate 1) 


Pointed (candidate 2) 


Search area (steradians) 


2.00x10"' 


2.35x10-5 


2.35x10-5 


log(Zl Jenkins) 


56351.1 


29692.1 


29687.3 


log(Zl.Evrard) 


56350.9 


29692.0 


29687.1 


log(Zo) 


56346.6 


29678.8 


29678.8 


Mv.Evrard 


0.11 


0.14 


0.14 


M^iJenkins 


0.29 


0.34 


0.34 




8.7 


7.9x10"* 


560 


^Jenkins 


26 


2.1 xio' 


1800 



to constrain cluster parameters such as P,rc,Mj 20Q ™d z. Using 
the Bayesian evidences we have calculated formal probabilities of 
detection taking into account point sources, receiver noise and the 
statistical properties of the primary CMB anisotropy. For the deeper 
component we find a formal probability of detection ratio of 7.9 x 
10'*:1 when assuming the Evrard et al. (2002) cluster number count 
and 2.1 X 10^:1 when assuming Jenkins et al. (2001) as the true 
prior. We derive a cluster mass 

• A Bayesian analysis using a phenomenological model of the 
gas distribution was also used to quantify the significance of the de- 
tection and again taking into account point sources, receiver noise 
and the statistical properties of the primary CMB anisotropy. For 
the deeper component we find ATq = — 295^[5/iK. 

• In our pointed follow-up observation the cluster system is de- 
tected with a high significance, with each map indicating that there 
is a 0.6mJy/beam peak decrement (8(TsA,pointed) towards the deeper 
component and an integrated decrement flux density (Ssz integrated) 
of ~ 1.2mJy/beam. The other component has a0.5mjy peak decre- 
ment and an integrated decrement of 0.7mJy. 

5/3 

• Using the approximation M^p oc Ssz.integrated we anticipate 
that the AMI blind cluster survey will detect clusters with Mj- 200 > 

2 X lO'^/i-g'Mg at4CTsA,sm-vey 
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Figure 9. ID and 2D marginal posterior distributions for the parameters in the SZ decrement model for candidate 1 (left) and candidate 2 (right). 
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